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Abstract 

We calculate the specific heats of quasi-particles of two-flavor QCD 
in its crystalline phases for low temperature. We show that for the 
different crystalline structures considered here there are gapless modes 
contributing linearly in temperature to the specific heat. We evaluate 
also the phonon contributions which are cubic in temperature. These 
features might be relevant for compact stars with an inner shell in a 
color superconducting crystalline phase. 



1 Introduction 



A number of theoretical studies have recently been devoted to QCD at low 
temperatures T and high densities. Besides the theoretical interest for the 
different QCD phases, the possibility of applications to compact stars, where 
dense quark matter might exist, has driven much of the recent interest. 

At high density and small T, quarks at the Fermi surface are expected to 
condense, giving rise to color superconductivity, see [TJ |2] and, for reviews, 

^On leave from the Dipartimento di Fisica, Universita' di Firenze, 1-50019 Firenze, 
Italia 
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At the largest densities QCD with three flavors should be in the so called 
CFL (color-flavor locked) phase. For decreasing density one expects pairing 
between quarks of non vanishing total momentum, resulting in crystalline gap 
behaviors [UElinilZI (for reviews see [E,^). Crystalline phases had already 
been discussed for electric superconductors (lUllIIl, so that one sometimes 
also refers to the crystalline phases as LOFF phases. A crystalline color 
superconductor shell might exist inside compact stars, between the external 
hadronic crust and the internal CFL core. 

In the present letter we discuss dispersion laws and speciflc heats of the 
quasi-particles present in the LOFF phase. These properties are relevant 
for the calculation of the thermal conductivity and neutrino emissivity and 
would affect the cooling of the compact stars. 

Thus far crystalline color superconductivity has been studied only in a 
two-flavor model. Different crystal structures have been computed, with the 
conclusion, based on a Ginzburg-Landau discussion, that the most favored 
structure at zero temperature is a cubic structure very close to a face centered 
cube (FCC). However at non zero temperature the situation is far from being 
clear (see the discussion in Ref. jHj), therefore we will discuss other crystalline 
structures as well. 

The results that one can easily predict for the speciflc heat in the LOFF 
phase with two flavors are a behavior for the phonons and a fi^T be- 
havior for the unpaired quarks. Our subsequent analysis (Sections 3 and 4) 
conflrms this expectation. For the paired quarks one would naively expect 
an exponential suppression exp(— A/T), which is the behavior of the regular 
2SC phase. However our analysis shows (see Sect. 3) that this is not the 
case, and for quark flavors that undergo LOFF pairing the speciflc heat has 
a /i^T dependence. The reason for this behavior is in the fermi quasi-particle 
dispersion law (see Sect. 2) that for the LOFF phase is gapless. As a con- 
clusion, in a neutron star, provided a LOFF phase is present, the speciflc 
heat and the thermal conductivity should be dominated by fermions, with 
paired and unpaired quarks on the same footing. A similar dominance is 
expected also for the LOFF phase with three flavors. In this case there are 
no unpaired quarks and therefore a a clear distinction should emerge in the 
thermal properties (speciflc heat, thermal conductivity) of the homogeneous 
case (CFL) and the inhomogeneous one (LOFF) because all the quarks in 
the former case would have the exponential suppression, while in the latter 
the T dependence is expected to be linear. 

The plan of the paper is as follows. In section 2 we discuss the disper- 
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sion laws of the fermionic quasi-particles for different crystal structures of 
the LOFF phase of high density QCD. We discuss the one-plane wave struc- 
ture, the strip, and the face centered cube. In section 3 we calculate the 
contribution of the fermionic quasi-particles to the specific heat for the three 
mentioned structures. In section 4 we use the effective phonon Lagrangian to 
calculate the phonon specific heats. Section 5 is devoted to the conclusions. 



2 Fermi quasi-particle dispersion law 

In this section we derive the Fermi quasi-particle dispersion law in the two- 
flavor QCD LOFF phase for a few different crystalline structures. We will 
work with an inhomogeneous condensate given by 

A(r) = A(r) 6,^3 6,, (1) 

where a, (3 are color indices and j = 1,2 are flavor indices. Notice that the 
gap term pairs together the color 1 and 2, whereas the fermions with color 3 
are unpaired. We work in presence of a difference in the chemical potential 
of the u and the d quarks. We define 

/i„ = /i + 5/i , /id = /i - . (2) 

As usual the quasi-particle dispersion law is obtained by looking at the zeros 
of the inverse propagator. For fermionic quasi-particles one has to solve the 
eigenvalue equation 

{S-')%x',=^. (3) 

where is the inverse propagator in the LOFF phase in the Nambu-Gorkov 
formalism and Xj are the Green eigenfunctions. Following the notations in 
[H] we write as follows: 



(4) 

where E is the quasi-particle energy and v is the Fermi velocity, that in QCD 
with massless quarks satisfies f = |v| = L Let us define 
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Performing the unitary transformation 

Qa ^ (^^is^ as ^-r/^' j , = (^e-'^" da) , (6) 

one measures the energy of each flavor from its Fermi energy. The resulting 
equations for F" and G° are independent of color and flavor indices, that 
therefore will be omitted below: 

(E + iv- V)G'-iA(r)F = 0, 

(E-iv V)F + iA(r)*G' = 0. (7) 

We solve these equations for three different crystalline structures, correspond- 
ing to different decompositions of (jT]) in plane waves: 

1. One plane wave; 

2. Two antipodal plane waves; 

3. The face centered cubic (FCC) structure, formed by eight plane waves 
with momenta pointing to the vertices of a cube. 

The first case is discussed for metals in |llj by Fulde and Ferrel (FF); for 
QCD is analyzed in |3]. The second case has been examined at T = by pH] : 
we will refer to it as the strip below. The last case is the preferred solution at 
T = in Ginzburg-Landau analysis of Ref. 7J. The reason why we examine 
these different structures is in the fact that the preferred structure at T 7^ 
can be different from that at T = 0, see e.g ^ and the discussion in 0. 

2.1 One plane wave 

For the FF condensate we take the direction of the Cooper pair total mo- 
mentum 2q along the axis. The gap 

A(r) = Ae^^^" (8) 

is therefore a complex number. If we take G = G'(r)e*^ '", F = F{r)e~^'^'^ , we 
get from Eqs.Q 

[E - qv, + w ■ V] G(r) = + iAF{r) , 

[E - qv, - zv ■ V] F(r) = - tAG{v) . (9) 
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These are the standard Gorkov equations for a uniform superconductor with 
energy E — qVz- The eigenfunctions are simple plane waves 

G{r) = Me'*^-"- F(r) = we'*^-"" , (10) 

and the quasi-particle spectrum is given by: 

E^ = qvz±^e + ^', (11) 

where ^ = k ■ v is the residual longitudinal momentum, i.e. the longitudinal 
momentum measured from the Fermi surface. Because of the transformation 
quasi-particle energies are computed from the corresponding Fermi en- 
ergies fiu,d- Eq. (fTT|l is the dispersion law of quasi-particle {E± > 0) or hole 
states {E± < 0). We notice that in this case, as in [T2], [IHl, [H] and [T3j , 
we are in presence of gapless superconductivity. 

An anisotropic dispersion law was also obtained in i4j by a different proce- 
dure (variational method). Their result was obtained for finite fi and reduces 
to (jllll if one considers only the leading order in the asymptotic /i — oo limit, 
which is the limit we consider in the present paper. 



2.2 Strip and FCC structures 

Both the strip and the cubic structure have real A(r). The solutions of Eqs. 
([7j) are Bloch functions 

G = n(r)e'''" , F = w(r)e'''" , (12) 

with u{r) and w{r) periodic functions and k in the first Brillouin zone. They 
satisfy 

[E - ^ + ■V]u{r) = +iA{r)w{r) , 

[E + ^- zv- V]w(r) = -iA{v)u{v) . (13) 

The corresponding quasi-particles are gapless [TU]. In fact, for = and 
^ = 0, the system (fT^ has two solutions. We have w± = ±u±, with u± 
solutions of 

V- Vu±(r) = ±A(r)M±(r), (14) 

and given by 

,,rf(r'-v) 



u±[r) = exp 



± / Afr', , 



(15) 
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Here the integration is over a path joining the origin to the point r. 

To find the dispersion law of quasi-particles for small values of ^ one uses 
degenerate perturbation theory and gets 



with 



A, 



1 





e 


A+A. 


/ (ir*exp 

/cell 





(/(r'- v) 



(16) 



(17) 



where Vc is the volume of a unit cell of the lattice. 

Let us now specialize to the case of the strip, i.e. the crystalline structure 
formed by two plane waves with wave vectors ±2q. The inhomogeneous gap 
is given by 

A(r) = A cos 2qz (18) 

and one gets 



At^ = yl(^) 



A 



(19) 



\ qv COS t) J 

where Io{z) is the modified Bessel function of the zeroth order. 

Let us now turn to the FCC structure. Summing the eight plane waves, 
the condensate can be put in the form 



A(r) = A cos 2qx cos 2qy cos 2qz 



(20) 



If we define 



B 



V ( sin 2q{x + y + z) sin 2g(x + y — z) sin 2q{x — y + z) 



+ 



sin 2q{—x + y + z) 



+ 



Vx+Vy- Vz 



-Vx + Vy + V^ 

one obtains for the cube 



^(fcc) 



i)7 dVe.J±^B] 

TX I J cell \ (l"^ ) 



(21) 



(22) 



where the integration is over the elementary cell of volume {T[/qY. 

Let us notice that whereas for the strip in Eq. ()19|1 one has a closed 
formula, for the FCC only numerical or approximate expressions can be given. 
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Since the parameter a; = A/gt> is expected to be small, expanding A^^'^' as a 
power series one obtains 

^/a^^%!)aF\^) = 1 + 0[x]\ (23) 

In conclusion for both the strip and the FCC the quasi-fermion spectrum is 
gapless. This result remains valid also for massive quarks, since the effect 
of the quark mass can be accounted for by reducing the quark velocity from 
the value v = 1 valid for the massless case jl6j. 



3 Specific heat of the Fermi quasi-particles 

The contribution of the Fermi quasi-particles to the specific heat per unit 
volume is 

dn f,^^dn{E,T) 



where, for the two flavor case (and |v| = 1), 

P=^, (25) 

while n{E, T) is the Fermi distribution function and the angular integration 
is over the directions of v. 

We should also remind that, in the case of two flavors, the two quarks 
of color 3 are unpaired and as such each contributes to the specific heat as 
follows: ^ 

Cv = ^T. (26) 

Let us stress however that quarks that do not participate in the standard 
spin condensate can nevertheless couple with spin 1, with much smaller 
gaps. Examples of these models can be found in [T7] . 

Let us now specialize Eq. pij) to the three crystalline structure under 
scrutiny. We limit the analysis to the small T range. 



3.1 One plane wave 

The dispersion law of quasi-particles is given by (fTTjl . Using one gets in 
the small temperature limit (T ^ A) and for A < q 

= '^f^ (q^-^ks) . (27) 



7 



The specific heat depends hnearly on temperature because the quasi-particle 
dispersion law (lllj) gives rise to gapless modes. There is also a contribution 
to the specific heat that comes from gapped modes, but this contribution is 
exponentially suppressed with the temperature. 

3.2 Strip and FCC structures 

From Eq. ()24|1 . using (fTT)|l that is valid for the strip and the FCC alike, one 
obtains, for T ^ A: 




This expression can be evaluated in closed form for the strip, when A± = A^^\ 
One gets 

c^:^ = ^ 1F2 (-1/2; 1/2, 1; {A/{qv)f) (quarks) (29) 

where 1F2 denotes the generalized hypergeometric function jTB]. Differ- 
ently from the analysis of [lOj, here v is not small and we can take z = 
A/qv near the second order phase transition. Since for small z one has 
iF2(— 1/2; 1/2, 1; z^) ^ 1 — z^, it is easily seen that the normal Fermi liquid 
result ()26|) is obtained for z = 0. On the other hand, at finite A, the specific 
heat turns out to be smaller. 

For the face centered cube, A^^^^ can be computed by evaluating the 
power expansion. One gets for this case 

cf ^) = ^ / ^./^p^ (quarks) . (30) 

Using the expansion ()23|) one can get an approximate formula giving a power 
series of A/qv. 

4 Specific heat of phonons 

Besides Fermi quasi-particles we consider also the massless Nambu-Goldstone 
bosons (NGB). Even though we expect a parametrically smaller 0{T^) con- 
tribution in this case, the contribution of the NGB might be relevant for 
future applications, e.g. for the calculation of thermal conductivity. For two 
fiavor QCD in the crystalline phase the only NGB are phonons ED! • 
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4.1 One or two plane waves 

Let us begin with the phonon Lagrangian for the plane wave 

1 



2(r-^y(V|i0r-<|Vx0|^j , (31) 

where is the phonon field, Vy = n ■ V, V_l = V — nV|| and n along the 
axis, i.e. the direction of q. In [1^] we found 

vf^ = cos^ , = ^ ^^^^ ^1 ' "^^2) 

with cos6'^ = 8112/1 = 0.833. Here 5/i2 is the 5/i value corresponding to 
the second order phase transition from the LOFF phase to the normal one. 
The dispersion law, relating the phonon quasi- momentum k and energy uo, 
therefore is 

a;(k) = ^vl{kl + kl)+vlkl . (33) 

Let be the total number of oscillatory modes. There is an oscillator for 
any quark pair and the total number of quark pairs is given by the available 
phase space, i.e. the pairing region. If V is the available volume, then 



V 2 J (27r)3' 

where g = 4 x 2 and we divide by two due to the fact that there are two 
quarks in the pair. The integration over the pairing region gives 

^-■^^"^'^(25). (35) 



V 2 {2n) 
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A factor 25 arises from the integration over the longitudinal residual momen- 
tum ^; 5 is the ultraviolet cutoff, of the order of fi while ( is the fraction of 
the phase space available for pairing in the LOFF phase. It is estimated of 
the order of A/q. 

The ratio N/V can be expressed in terms of a cutoff frequency by a 
procedure analogous to the introduction of the Debye frequency for ordinary 
crystals. We write 

- = f{u)du: (36) 
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where /(cu) is given by 



Substituting in (jMj) we see that cuo is of the order of /^(A/g)^/'^ therefore 
large (excluding the second order phase transition region). More precisely 
we get the formula 

N gujl 



V 127r2t;2' 



(38) 



These results also hold in the case of the strip, i.e. two antipodal plane 
waves. As a matter of fact one can prove, following the same procedure as in 
that the effective Lagrangian for the phonon field is still given by 

The specific heat per unit volume at small temperatures {uo ^ T), is 
given by {ks = h = 1): 

—V (39) 



5V \ujdJ 
and therefore 

c(FF) = = (phonons) . (40) 

I5vj_v\\ 

This result holds for the two considered structures (one or two plane waves). 
4.2 Face-centered-cube 

For the FCC crystal structure the phonon Lagrangian is given by [21] 

^ 4=1,2,3 j=l,2,3 ^ j=l,2,3 i<i=l,2,3 

(41) 

In |20|| we found 

a = 1/12 6 = 0, c= (3cos2 0g-l)/12. (42) 
Eq. (j36|) is still valid with 

/M = 5i:f/(g^(--.W). (43) 
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since u!r{k) = Vr{'n)k we get 

f(uj) = — ^ > / -— — = — ^ K . (44) 

Here the velocities v"^ = vf.{n) are the eigenvalues of the matrix 

I a + bnf cnin2 criin^ \ 

cnin2 a +hn\ cn2n3 , (45) 
\ criin^, cn2n3 a +hn\j 

with rii = sin ^ cosy?, n2 = sin^siny^, = cos 6. 

To get the Debye frequency and the specific heat we use the numerical 
result 

K=^i:fJ^''3^3ylO' (46) 



corresponding to the values (j42j) of the parameters. For the specific heat we 
get, similarly to Eq. pn|) . the result 

cf^)= iTT^' (phonons) . (47) 

One can note that numerically c^^^^ > 



5 Conclusions and Outlook 

We have considered three different crystalline structures for high density 
QCD in the LOFF phase: one plane wave, the strip, and the favorite face- 
centered-cubic strucure. Gap less fermionic quasi-particles are found to exist 
even in the presence of quark mass terms. These gapless fermions provide for 
the dominant contributions to the thermal properties for low temperatures, 
such as of possible interest for compact stars. Typically they contribute to 
specific heats with terms linear in T which are dominant at low T. The 
massless phonons contribute instead as while gapped particles in the ho- 
mogeneous case (2SC) have an exponential suppression at low T. The pre- 
vious calculations could be applied to the evaluation of transport properties 
of compact stars. At the present a reliable model for the LOFF phase in 
compact stars is still lacking, since the study of the inhomogeneous super- 
conducting phase of QCD with three flavors has not yet been performed. 
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Nevertheless one can make some conjecture by imagining the existence of 
the LOFF phase inside the star, but outside the inner core, where one might 
suppose the existence of quark matter in the CFL phase. In this analysis the 
thermal properties of the homogeneous and inhomogeneous phases would be 
different since for the homogeneous case the fermion specific heats would be 
exponentially suppressed, while for the inhomogeneous case the vanishing 
with the temperature would be linear. An application of previous results 
might be the evaluation of the thermal conductivity in the LOFF phase for 
the existing model of 2-fiavor QCD. We plan to come to this subject in the 
future. 
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